PReBeaM for Planck: A Polarized Regularized Beam 
Deconvolution Map-Making Method 
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ABSTRACT 

We describe a maximum likelihood regularized beam deconvolution map-making algorithm 
for data from high resolution, polarization sensitive instruments, such as the Planck data set. 
The resulting algorithm, which we call PReBeaM, is pixel-free and solves for the map directly 
in spherical harmonic space, avoiding pixelization artifacts. While Fourier methods like ours are 
expected to work best when applied to smooth, large-scale asymmetric beam systematics (such 
as far-side lobe effects) we show that our ?7i-truncated spherical harmonic representation of the 
beam results in negligible reconstruction error - even for m as small as 4 for a polarized elliptically 
asymmetric beam. We describe a hybrid OpenMP/MPI parallelization scheme which allows us to 
store and manipulate the time-ordered data from instruments with arbitrary scanning strategy. 
Finally, we apply our technique to noisy data and show that it succeeds in removing visible power 
spectrum artifacts without generating excess noise on small scales. 

Subject headings: cosmic microwave background - cosmology: observations - methods: data analysis 



1. Introduction 

One of the most exciting prospects for the up- 
coming Planck satellite is its capability to mea- 
sure the polarization anisotropics of the CMB over 
the entire sky in nine frequency channels. The po- 
tential rewards from these measurements are many 
and include tighter constraints on cosmological pa- 
rameters, determination of the reionization history 
of the universe, and detection of signatures left by 
primordi a.1 gravitational waves generated during 
inflation (JThe Planck Collaborationll2005[) . 

Measurement of the CMB polarization signal 
presents a great experimental challenge as it is an 
order of magnitude smaller than the temperature 
signal and is especially susceptible to distortions 
due to optical systematics and foreground contam- 
inants. Indeed, if left untreated, leakage from the 
much stronger temperature signal will contami- 



nate the polarization maps. Maps and spectra will 
also suffer from leakage from E-mode polarization 
to B-mode polarization, jeopardizing the potential 
detection of inflationary B-modes. At the resolu- 
tion and sensitivity of the next generation of ex- 
periments, including the Planck mission, studies 
of primordial non-Gaussianity may also be sensi- 
tive to beam-induced systematics. In this paper, 
we present a novel technique for both assessing 
and removing systematic effects due to beams in 
temperature and polarization maps. 

The Planck satellite is designed to extract 
essentially all of the information in the pri- 
mordial temperature anisotropics and to mea- 
sure the polarization anisotropics to high accu- 
racy for 2 < £ < 2500. This wiU be achieved 
by measuring the full-sky signal to an angu- 
lar resolution of 5', to a sensitivity of AT/T 
^ 2xl0~^, and over a frequency range of 30-857 
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GHz (|The Planck Col laboration 2005). The sci- 
entific performance of Planck depends, in part, 
on the behavior of systematic effects which may 
distort the signal. 

A primary objective of Planck is to produce 
all-sky CMB maps at each frequency. The pro- 



cess by which the satelhte's time-ordered data 
(TOD) is wrapped back on to the sphere to cre- 
ate an image is known as map-making. The map- 
making process becomes difficult due to a num- 
ber of chaUenges: distortions in the beam, fore- 
ground contamination through far-side lobes, size 
of the data, and correlated noise effects. It is of 
critical importance to fully characterize the beam, 
and use this information during map-making to 
deconvolve beam effects. We have previously de- 
scribed a powerful map-making algorithm which 
implements the beam deconvolution technique for 
temp erature measurements (JArmitage fc Wandeltl 
20041 ). In this paper, we will extend that de- 
scription to include polarization measurements. 
We refer to this new technique as PReBeaM: 
Polarized Regularized Beam deconvolution Map- 
making. While we focus on reconstructing the 
map with a uniform effective beam and realize 
corrections to the p ower spectrum as a co nse- 
quence, other work bv lSouradeep et al.l (|2006l ) and 
Mitra et al.l (|2007l) has focused on deriving correc- 
tions to the power spectrum due to asymmetric 
(non-circular) beam effects. 

Within the Planck collaboration, the CTP 
working group has developed five map-making 
methods and compared their results using the sim- 
ulated 30 GH z data in what is kn own as the 
Trieste paper (JAshdown et al.l l2008h . The Tri- 
este paper assessed the impact of beam asymme- 
tries on the Planck spectra without attempt- 
ing to treat the problem of beam asymmetry at 
the map-making level (an angular power spectrum 
correction method was developed based on sim- 
plifying assumptions). In addition to PReBeaM, 
another deconvolution map-makin g technique for 



Planck has been established by [Harrison et al. 
( 20081 ). Both methods allow for arbitrary beam 



shapes and in both cases the asymmetry of the 
beam is parametrized by an asymmetry param- 
eter rrirnax which can vary between and tmax- 
Our method scales computationally as a function 
of rrimax] this is advantageous when large gains in 
accuracy can be achieved with small increases in 
iTT-max- In contrast, the Harrison method incurs a 
fixed computational expense for arbitrarily large 
iTimax- The Harrison method takes advantage of 
the Planck scanning strategy to condense the full 
TOD into phase-binned rings, thereby achieving a 
significant reduction in processing time. 



A complete characterization of the beam in- 
cludes both the main beam and the far-side lobes. 
Sidelobes are located as far away as 90° from the 
main focal plane beam, and therefore require a 
large rrimax par ameter for a complete harmonic 
description. In lArmitage &: Wandeltl (|2004l ) we 
demonstrated the full potential of our method us- 
ing far-side lobes and maps with foreground sig- 
nals. Here, we show the usefulness of PReBeaM 
for deconvolving main-beam distortions. In fact, 
we find that it makes sense to use PReBeaM for 
main beam effects since only a small rUmax param- 
eter is needed to capture the azimuthal structure 
of the main beam. In this way, we profit from 
the computational advantage of our method in the 
case of small rrimax, allowing for the unified treat- 
ment of main beam and side lobe effects. 

In ^21 we describe the deconvolution map- 
making algorithm for PReBeaM. The simulated 
data and beams are detailed in fJ31 We present 
results in 21 showing the effectiveness of PRe- 
BeaM in removing systematic effects due to beam 
asymmetry and we discuss computational consid- 
erations. We finish with our conclusions from this 
study in SJS) 

2. PReBeaM Method 

First we review the standard set-up to the map- 
making problem for a solution of the least-squares 
(or maximum-likelihood) type. 

The TOD generated by a detector is effectively 
a convolution of the true CMB sky with a beam 
function. If we consider the sky as a pixelized vec- 
tor, it will have length ripixei x ripoi where ripoi = 3 
for the I (total intensity), Q, and U Stokes com- 
ponents. The riTOD-length TOD vector d is the 
result of a matrix multiplication of the observation 
matrix A with the sky s 



As = d. 



(1) 



In our implementation of the maximum- likelihood 
solution, we refer to A as the convolution opera- 
tor. A encodes information about both the scan- 
ning strategy and the optics of the scanning in- 
strument. The least-squares estimate of the true 
sky, s, is given by the normal equation 

A^As = A'^d (2) 

where A""" is the transpose convolution operator. 
Equation Q is exact if the noise is stationary 



and uncorrelated in the time-ordered domain. The 
generahzation to non-white noise is as follows 



A^N-^As = A'^N-M 



(3) 



where N is a noise covariance matrix. In this work 
we consider CMB only and CMB plus white noise. 
We modify the normal equation by introduc- 
ing a regularization technique in order cope with 
the ill-conditioned nature of the coefhcient matrix 
A"^ A. We split off the ill-conditioned part of A by 
factoring it into two parts: A — BG. The factor G 
is what we refer to as the regularizer in PReBeaM. 
In general, the regularizer can be any target beam; 
a natural choice would be the angle-averaged de- 
tector beam. In our study, we choose G to be a 
Gaussian smoothing matrix, defined in harmonic 
space as 



G\ 



G 



G,C 



exp 



exp 



-(j^(.{t+l) 



-a^{(.{i+l)-A) 



(4) 



where a = FWIIM/-\/81n2. The superscripts G 
and G refer to the gradient and curl components 
in the typical linear polarization decomposition. 
Our modified normal equation becomes 

B^^Bx = B'^d (5) 

where we are solving for x = Gs. In this way, we 
do not attempt to reconstruct the sky at a higher 
resolution than that of the instrument. 

In a standard pixel-based solution of equation 
([2|), in which one assumes that the observing 
beam is spherically symmetric, A is a sparsely- 
filled pointing matrix. For polarization measure- 
ments, each row of A contains only three non- 
zero elements. The deconvolution map-making 
approach does not assume spherically symmet- 
ric beams, instead allowing for arbitrary beam 
shapes. We achieve this added complexity pri- 
marily by solving the normal equation in spherical 
harmonic space in order to make use of fast and 
exact algorithms for the convolution and trans- 
pose co nvolution of two arbitrary functions on the 



spher e (|Wandelt fc Gorskil 120011 IChallinor et al. 
20001 ). These algorithms are described in abbre- 



viated form below. A secondary advantage of op- 
erating entirely in harmonic space is that artifacts 
due to pixelization (such as uneven sampling of 
the pixel) are completely avoided. 



2.1. Fast all-sky convolution for polari- 
ametry measurements 

For a full presentation of the formalism for con- 
volution of an instrume nt beam with a s ky signal, 
the reader is referred to lChallinor et al.l ( 20001 ). 

In compact spherical harmonic basis, equation 
^ is written as 



-^L' M' mm'm" -'- mm' m" 



^L' M'mm'm"^mm'm" LMSLM 

(6) 
where s^m is the spherical harmonic representa- 
tion of the sky and Tmm'm" is defined as the re- 
sult of a convolution of a band-limited function 
b wit h the sky s. The Planck Level-S soft- 
ware ( Reinecke et al.ll2006l ) nomenclature refers to 
Tmm'm" as a ring set. This is written in harmonic 
space as 

J-mm'm" — 2^\r,^em"eM' + ^em'^lM' 



+ S£mbeM')dmM (^B)c^Af M' W 



(7) 



where {9e,0) are fixed parameters which define 
the scanning geometry. 

In equation ([7]), d^mAii^E) and d^MM'i^) ^'''^ ^^" 
lated to the Wigner D-matrices by 



(8) 



Analogously, the transpose convolution in har- 
monic space is given by 



Vim 



m'm' 



dmm'{(^E)dm'm"{(^)be^"T, 



mm'm" 



(9) 



where P = I,G, G. 

2.2. PReBeaM Implementation 

Now we outline the algorithmic steps taken to 
make a map from a TOD vector by PReBeaM. 

First we construct the right-hand side of equa- 
tion ([6]) in two steps: converting TOD to a 



T, 



mm'm' 



array and applying A . T, 



mm'm" 



structed by transpose interpolating the TOD vec- 
tor d. The transpose interpolation of the TOD 
vector onto the Tmm'm" grid is akin to a binning 
step, where each element of the TOD is mapped, 
via interpolation weights, to several elements of 
the Tmm'm" cubc according to the orientation 
and position of that data point in the scanning- 
strategy. The interpolation scheme is described in 



greater detail in ^2.31 Next, we transpose convolve 
the beam coefficients bim with Tmm'm" according 
to equation ^. 

Once the right-hand side has been computed, 
we use the conjugate gradient iterative method to 
solve equation ([2])- With each iteration, the coef- 
ficient matrix A'^A, is applied using the following 
procedure: 

1. Apply the convolution operator. A, to 
project the sky a^m on to the convolution 

gnu J-mm'm" 

2. Inverse Fourier transform over the first two 
indices of Tmm'm" to get r$2,e,m" (we omit 
the transform over m" as it is incorporated 
in the interpolation scheme) 

3. Forward interpolate from T^2,Q,rn" to a 
TOD vector 

4. Transpose interpolate from the TOD vector 
to a new ring set T'4^_e,m" 

5. Fourier transform over the first two indices 

6. Apply the transpose convolution operator, 
A"^, to project the ring set T'^^,^„ back 
into a new sky &£„ vector 

2.3. Polynomial Interpolation and Zero- 
Padding 

PReBeaM uses the same polynomial interpola- 
tion as implemented in the Level- S software used 
to generate the sirnulati on TODs and as described 



Reinecke et alj (|2006[ ). The objective of forward 



interpolation is to construct a TOD element at a 
particular co-latitude, longitude and beam orien- 
tation using several elements of the ring set T and 
their corresponding weights. Transpose interpola- 
tion operates in exactly the opposite manner as 
the forward interpolation: distributing a single el- 
ement in the TOD to multiple elements of the ring 
set. This is done using the same weights calculated 
for the forward interpolation. The entire opera- 
tion of interpolation and transpose interpolation 
from ring set to TOD and back again is depicted 
in figure [1] 

PReBeaM also includes the option to zero-pad 
during the FFT and inverse FFT steps. This 
means that the working array (either Tmm'm" or 
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Fig. 1. — Forward interpolation from ring set to 
TOD element and transpose interpolation from 
TOD element to ring set. 

T^2,0-m") is enlarged and padded with zeroes out 
to imax.pad > f^max- This has the cffect of decreas- 
ing the sampling interval. We found that the com- 
bined effects of small-order polynomial interpola- 
tion (order 1 or 3) and zero-padding of 2 x £max 
or 4 X (max dramatically reduced the residuals in 
our maps. 

2.4. Parallelization Description 

PReBeaM employs a hierarchical paralleliza- 
tion scheme using both shared-memory (OpenMP) 
and distributed-memory (MPI) types of paral- 
lelization. The map-making was performed on the 
NERSC computer Bassi. Bassi processors are dis- 
tributed among compute nodes, with 8 processors 
per node. OpenMP tasks occur within a node and 
MPI tasks occur between nodes. 

We show a diagram of our hybrid paralleliza- 
tion scheme in Figure [D The full TOD and point- 
ings are divided equally between the nodes for in- 
put and storage of pointings. Within an iteration 
loop, four head nodes are designated to perform 
the convolutions, while the remaining active nodes 
are dedicated to the interpolation routines. Each 
of these four nodes performs the convolution of 
the sky with one of the four detectors. The re- 
sulting arrays are then distributed to all nodes for 
interpolation over the segment of data stored there 
and then gathered back onto the designated nodes 
for transpose convolution. Finally, the agm are 
summed, using MPI task mpi_reduce, into a sin- 
gle aim on a single node; this is the new estimate 
for the sky vector. 

This particular scheme was devised so that the 
four distinct convolutions that must occur (one 



sky with four different beams) can take place si- 
multaneously, while the pointings are distributed 
among as many nodes as possible for maximum 
speed in interpolation. Both convolution and in- 
terpolation and their transpose operations make 
use of all processors on a node by using OpenMP 
directives. 




convolution 
with beam 1 



interpolation and 
transpose 
interpolation on 
TOD segment 



convolution 
with beam 2 




interpolation and 
transpose 
interpolation on 
TOD segment 




Fig. 2. — Depiction of hybrid parallelization 
scheme used in PReBeaM. Rectangles represent 
work done on a node, ellipses represent data prod- 
ucts, and arrows represent transfer of data. The 
work done within a node (convolution, interpola- 
tion and their transpose operations) is parallelized 
using OpenMP. This shows a slice of two head 
nodes though the algorithm may operate on many 
more nodes. 



3. Simulations and Beams 

The simulated Planck data on which PRe- 
BeaM was run was generated by the Planck CTP 
working group for the study of the performance 
and accuracy of five map- making codes summa - 
rized in the Trieste paper ( Ashdown et al.ll200a l. 
Planck will spin at a rate of approximately one 
rpm, with an angle between the spin axis and the 
optical axis of ^ 85°. We used the cycloidal scan 
strategy in which the spin axis follows a circular 
path with a period of six months, and the angle 
between the spin axis and the anti-Sun direction 
is 7.5°. TODs were generated for 366 days for the 
four 30 GHz Low Frenquency Instrument (LFI) 
detectors. At a sampling frequency of 32.5 Hz, 
this corresponds to 1.028x10^ samples per detec- 
tor, for a total of over 65 Gb of data and point- 
ings. The simulated data also included the effects 
of variable spin velocity and nutation (the option 
to include the effects of a finite sampling period 
was not included). 

The data was simulated with elliptical having 
with a geometric mean FWHM of 32. '1865 and el- 
lipticity (maximum FWHM divided by minimum 
FWHM) of 1.3562 and 1.3929 for each pair of 
horns. The widths and orientations of the beams 
were different; this was referred to as beam mis- 
match in the Trieste paper. In spherical harmonic 
space, the simulation beams were described up to 
a beam rrimax of 14. The same beams were used 
in PReBeaM to solve for the map, although we 
allowed the beam asymmetry parameter rrimax to 
vary. 

4. Results and Discussion 



For this paper, we make temperature and po- 
larization maps from simulated one-year observa- 
tions of the four 30 GHz detectors of the Planck 
LFI. We examine two cases: CMB signal only and 
CMB plus uncorrelated (white) noise. Foreground 
signals and correlated noise properties will be ex- 



amined in a future paper. The 30 GHz data was 
an optimal choice for this analysis because the low 
sampling rate and resolution minimize the data 
volume, while the large beam ellipticity allows us 
to demonstrate the full potential of our beam de- 
convolution technique. 

PReBeaM operates entirely in harmonic space, 
solving for and producing as output aim- For vi- 
sualization purposes, maps were ma de from aim^s 
out to £max 512 and at the Healpix (JGorski et al. 



20051 ) resolution of nside 512 {~ 7" pixel size). 
Most of the results presented in this paper were 
attained with an FFT zero-padding of factor four, 
an interpolation order of 3, and an asymmetry pa- 
rameter of TTimax =4 (wc notc whcrc the param- 
eters differ from this). To compare with the in- 
put signal, a reference map representing the true 
sky was created by smoothing the input aim by 
a Gaussian beam of FWHM = 32. '1865. Sim- 
ilarly, our regularizer G (in equation ([¥])) was 
set to have a FWHM of 32. '1865 to match this 
smoothing. As noted in fJS] the same data we 
use here has been processed by five map-making 
codes in lAshdown et al. (2008). We have chosen 
to compare our results with the analogous results 
from Springtide, one of the codes in this study. 
Springtide was chosen, out of the five codes, be- 
cause it is the map-making code installed in and 
used by the Planck Data Processing Centers for 
the HFI and LFI instrument. It is sufficient to 
compare with Springtide only as no significant 
differences in accuracy were found between codes 
(with similar baselines and in the absence of noise) 



(jAshdown et al.l 120081 ) . In the absence of noise, 



Springtide is algorithmically akin to a straight- 
forward binning of the TOD into a sky pixel map. 
Because we are using Springtide to represent all 
non-beam-deconvolution methods we will refer to 
the Springtide maps as the binned maps. 

We begin by examining the spectra in the 
binned map, PReBeaM map and the smoothed 
input map shown in Figure [3l The effect of the 
beam mismatch is clearly seen where the peaks 
and valleys of the binned map spectra have been 
shifted towards higher multipoles. The detectors 
measure different Stokes I which translates to ar- 
tifacts in the polarization map. Deconvolution 
suppresses leakage from temperature to polariza- 
tion as evidenced by the PReBeaM spectra which 
overlays the input spectra. This shift is expected 



to remain apparent in the TE spectra of non- 
beam-deconvolved maps even in the presence of 
noise because of larger temperature signal and the 
temperature-to-polarization cross-coupling. 

The fractional difference in the angular power 
spectrum (defined as {Ce^^^ — C£i„)/C£i„) of the 
input and output maps is shown in Figure [H We 
show the fractional difference spectra for the TT, 
EE, and cross-correlation TE signals, omitting the 
BB spectra since Cf^ is zero in the simulation of 
the CMB map . The results for PReBeaM are 
shown at three intervals: the 25th, 50th and 75th 
iterations. This shows the behavior of the power 
spectra as PReBeaM converges on the solution. 
The beam mismatch effect is also seen in Figure 
m where the fractional difference in the PReBeaM 
spectra lie closer to zero than the binned map 
spectra over the full range of multipole moments 
for EE and TE. 

As described earlier, PReBeaM allows for vari- 
ation in the asymmetry parameter rrimax- We ex- 
amined the performance of PReBeaM as a func- 
tion of mmaxi setting it to 2, 4 and 6. A remark- 
able improvement in the power spectra was found 
by increasing mmax from 2 to 4, while an increase 
from 4 to 6 only resulted in marginal improve- 
ments. This effect is best seen in the BB power 
spectra as shown in Figure El Thus, while the in- 
put TOD was simulated with a beam having an 
"mmax cut-off of 14, PRcBcaM operates optimally 
at an mmax of just 4, thereby allowing us to cap- 
italize on the computational property that PRe- 
BeaM scales as rumax- 

We define a quantity called n^- 



E 

e'=2 



\ACe'\ 

^ Planckp/ 



(10) 



which we use to quantify the maximum, or worst- 
case bias beam systematics could induce in a cos- 
mological parameter that happened to be degener- 
ate with that parameter. The quantity upianck is 
the expected one-sigma errors for the LFI 30 GHz 
channel, computed as the diagonal elements of the 
covariance matrix for the simulated input spectra, 
assuming a sky fraction of 0.65. The n 



sigyna 



values 

are plotted in Figure [5] and show that PReBeaM 
reduces the worst case bias due to untreated beam 
systematics from tens of sigma to much less than 
one sigma over the entire H. range. 
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Fig. 3. — EE, TE, and BB spectra of smoothed input map (black curve), binned map (blue curve) and 
PReBeaM (red curve). The EE and TE spectra show the effect of temperature-to-polarization cross-coupling 
seen in the binned map spectra as shifts in the peaks and valleys and absent from the PReBeaM spectra. 
The input BB spectra is absent from the BB plot since the input B-modes were zero. TT spectra are omitted 
since differences in the three spectra are not apparent in this representation. 
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Fig. 4. — Fractional difference in power spectrum for PReBeaM (red, blue, and cyan curves) and the binned 
map (black curve) for TT, EE, and TE. Spectra for PReBeaM are shown as a function of number of iterations 
to demonstrate convergence. 
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Fig. 5. — BB power spectra as a function of asymmetry parameter rrimax for nimax = 2 (blue curve) , 4 
(cyan curve), and 6 (red curve). The input BB spectra was zero so the smallest output BB spectra is most 
desirable. In this run, the PReBeaM input parameters interpolation order and zero-padding were set to 1 
and 2, respectively. 
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Fig. 6. — Worst case bias in estimation of cosmological parameters due to errors in the power spectra of 
PReBeaM (dashed curve) and due to the errors in the power spectra of the binned map (solid curve). 



We examine the resulting temperature and po- 
larization (Q and U) maps. The output map 
for both PReBeaM and the binned map was sub- 
tracted from the smoothed input map at the same 
resolution to make the residual maps shown in Fig- 
ure [71 PReBeaM residuals were plotted on the 
same color scale as the binned map, showing that 
PReBeaM attained smaller residuals for both tem- 
perature and polarization. 

As a final test, we run PReBeaM on TOD con- 
taining CMB signal and white noise and compare 
with the smoothed input CMB spectrum and the 
analogous results from Springtide (in this case we 
refer to Springtide directly since this is not sim- 
ply a binned map). The level of the uncorrelated 
noise is specifed in the detector database and has 
a nominal stan dard deviation per sa mple time of 
a = l350fj.K (JAshdown et al.ll2008l ). PReBeaM 
achieves a noticeably superior fit to the input spec- 
trum compared with Springtide from £ ^ 150 to 
~ 250. Assessing the relative performance of PRe- 
BeaM and Springtide in more detail would require 
performing Monte Carlo averages. We focus on the 
TE spectrum since the improvement is visible even 
for a single simulation. For the other spectra PRe- 
BeaM performs as least as well as Springtide but 
the detailed difference are more difficult to assess 
without a Monte Carlo study. 

4.1. Computational Considerations 

The computational costs and advantages of our 
method should be noted. To perform a convo- 



requires 0{i^ 



max"^rnax 



) for the 



lution up to i 
general case. Since rrimax is bounded by imax , the 
cost never scales worse than ©(^^ax) ^^'^ is only 
^i^max) foi' the symmetric beam case. By com- 
parison, a brute force computation in pixel space 
would require 0{i^g^^). In this study, data was 
simulated with beams having an asymmetry pa- 
rameter of nimax = 14, but maps were made using 
a cut-off value of rrimax = 4 in PReBeaM. We have 
demonstrated that computational cost can be con- 
served while still achieving the benefits of beam 
deconvolution 

It was found that an increase in the zero- 
padding factor from two to four produced superior 
results over an increase in the interpolation order 
from one to three. An optimal run of PReBeaM 
will therefore include the largest zero-padding pos- 
sible given machine memory constraints in con- 



junction with a polynomial interpolation of or- 
der one or three. This is advantageous since 
the time spent in an FFT is nearly negligible 
and affected only minimally with an increase in 
zero-padding. In contrast, time for interpolation 
scales as interpolation-order-squared and as this 
is a TOD-handling step, it dominates over any 
cost incurred by convolutions. In the case of the 
results shown here, interpolation steps consume 
more than 90% of the wall-clock time per itera- 
tion. 

The results produced here were generated using 
12 nodes on NERSC computer Bassi (making use 
of all 8 processors per node) and was complete in 
about 29 wall-clock hours, for a total of 2797-CPU 
hours. The maximum task memory was 20 GB on 
a single node. 

5. Conclusion 

We have found that PReBeaM has outper- 
formed the standard binned noiseless map using 
two measures: spectra and residual maps. We 
examined the fractional differences in the spec- 
tra and found markedly smaller differences in the 
PReBeaM spectra versus the binned map spectra 
across a range of multipole moments. We find that 
map-making codes which do not deconvolve beam 
asymmetries lead to significant systematics in the 
polarization power spectra measurements. The 
temperature-to-polarization cross-coupling due to 
beam asymmetries is manifested as shifts in the 
peaks and valleys of the spectra. These shifts are 
absent from the PReBeaM spectra. We translated 
the errors found in the power spectra to an esti- 
mate of the statistical significance of the errors in a 
parameter estimation resulting from these spectra, 
which we call n^- This analysis showed that the 
worst case parameter bias due to beam-induced 
power spectrum systematics could be tens of sigma 
while PReBeaM reduces the risk of parameter bias 
due to beam systematics to much less than 1 sigma 
We also found the I, Q, and U component resid- 
ual maps to be smaller for PReBeaM than for the 
binned map, implying smaller map-making errors. 

We have presented here the first results from 
PReBeaM for a straightforward test cases of CMB 
only and CMB plus white noise, and including 
only the effects of beams in the main focal plane. 
However, there is great potential for using PRe- 
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PReBeaM Residual T 
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PReBeaM Residual Q 
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PReReaM Residual U 



Springtide Residual U 





Fig. 7. — The residuals between the input reference sky and PReBeaM output (left column) and the residual 
between the input reference sky and the binned map (right column) for Temperature (T), and the Stokes Q, 
and U parameters. 
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Fig. 8. — TE spectrum of CMB and white noise for Springtide (blu e curve) and PRcBeaM (red curve). The 
smoothed input map is shown in black. Following the example in I Ashdown et al., (,2008l) . we reduce i to i 
variation by filtering the spectra by a sliding average {A£ = 20). In this run, the PReBeaM input parameters 
interpolation order and zero-padding were set to 1 and 4, respectively. While PReBeaM performs at least as 
well as Springtide in the TT, EE, and BB spectra, we omit these spectra since the detailed differences are 
difficult to assess without an in-depth Monte Carlo study. 



BeaM to remove or assess systematics due to 
the combination of foregrounds and beam side 
lobes. Systematics introduced by side lobes will 
appear on the largest scales, potentially imped- 
ing the detection of primordial B-modes on the 
scales where they are most likely to be measured. 
We hav e already shown for tempera ture measure- 
ments (JArmitage &: Wandeltl l2004h that our de- 
convolution technique can be used to remove ef- 
fects due to side lobes. Future work will examine 
the noise properties of PReBeaM maps and will in- 
clude foregrounds from extragalactic sources and 
diffuse Galactic emission. 
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